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The effect of meson retardation in A^Ai'-interaction and exchange currents on deuteron electrodis- 
integration is studied in a coupled channel approach including NN-, NA- and yrd-channels. It is 
shown that the influence of retardation depends on the energy regime: Whereas below 7r-threshold 
calculations with static and retarded operators yield almost identical results, they differ significantly 
in the A-region. Especially, the longitudinal and the longitudinal-transverse interference structure 
functions are strongly affected. 
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I. INTRODUCTION 



o 
o 

(N 

Recently we have constructed a realistic iVA^-interaction model and corresponding electromagnetic (e.m.) exchange 

■ currents based on a coupled channel approach with meson, nucleon and A-degrees of freedom and have applied it 
to A^A^-scattering and deuteron photodisintegration Q, 0, 13 ■ It contains complete meson retardation in both the 
A'^iV-interaction as well as in the e.m. pionic meson exchange currents (tt-MEC). In addition off-shell contributions to 

K>" I the e.m. one-body current were considered in subsequent work '4'|. As was shown in the influence of retardation 

■ is large above pion threshold and only its incorporation leads to a satisfactory theoretical description of deuteron 
\^ photodisintegration in the A-region. On the other hand, e.m. off-shell effects turned out to be quite small ,43 . 

, Due to the possibility of an independent variation of energy and momentum transfer in the space-like region and 
' the existence of a longitudinal polarization of the exchanged virtual photon, it is expected that electrodisintegration 
■<!;;j- , of the deuteron will provide additional insights into the hadronic properties of the iV TV-system. For example, it 
' allows the study of the A'^A'^-interaction in the short-range regime independent of the excitation energy in contrast 
to photodisintegration. Moreover, in quasifree kinematics this reaction is an important tool for the exctraction of 
-j-J ' nucleon properties like, e.g. the electric form factor of the neutron For this reason, we have extended our 

1— H ' approach to electrodisintegration. In the next section, a brief survey of our model is given. Then, in Sect. UTTl results 
^ , are presented and compared to experimental data for some selected kinematics and, finally, some conclusions are 
^ drawn. 

> : 

; II. THEORETICAL FRAMEWORK 

H ; 

5t 1 The model Hilbert space consists of three basic configurations, i.e. two nucleons (NN), one nucleon and one delta 
(A'^A), and two nucleons and one meson (e.g. NNn). The basic hadronic interactions are determined by baryon- 
baryon-meson vertices, where we include as mesons tt, p, a, S, to, and rj. Inserting the vertices into the appropriate 
Lippmann-Schwinger equation, one obtains after some straightforward algebra [21 the desired retarded one-boson- 
exchange operators describing the interactions A^A^ ^ A^A^, A^A^ ^ A^ A, and A^A ^ NA. In our explicit realization, 
we use for the parametrization of the retarded A'^A'^- interaction the Elster potential j^J, which includes in addition 
one-pion loop diagrams as nucleon self energy contributions in order to fulfil unitarity above pion threshold. For this 
reason, one has to distinguish between bare and physical nucleons. For details we refer to [3- For the interactions 
A^A^ <-> A^A and NA ^ NA, we take besides retarded pion exchange in addition only static p-exchange into account. 
Moreover, the interaction of two nucleons in the channel with deuteron quantum numbers in the presence of a 
spectator pion, called 7rd-channel, is also considered. By a suitable box renormalization 0, an approximate phase 
equivalence between the Elster potential, which does not include A-d.o.f., and our coupled channel approach below 
pion threshold is obtained. The hadronic A-parameters are fitted to the P33-channel of pion-nucleon scattering and the 
^£)2-channel in A^A^-scattering. For the sake of comparison with the conventional approach with static interactions. 
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we have constructed also a static version of our model interaction based on the Bonn OBEPR-potential 1^ where in 
accordance with the usual treatment the pion-nucleon loop diagrams have been neglected. For the results presented 
below, we have used two of the various model versions discussed in 0, namely "CC(ret,7r,p,0)" and "CC(stat,7r,p,0)" 
for the retarded and static calculations, respectively. 

The basic e.m. interaction of the model is described in detail in 1^ Q| and consists of baryon and meson one-body 
currents as well as Kroll-Rudermann (contact) and vertex contributions. These currents are, together with the ttNN- 
vertex, the basic building blocks of the corresponding effective current operators for two baryons consisting of one- and 
two-body contributions. The matrix element of the effective nucleon one-body operator between plane- wave states of 
two bare nucleons N is given in J] for real photons. Its extension to virtual photons with four momentum ~ (w, q) 
reads as follows 

1^1 ^, i^7V[i]aj/ ^^|^ - \ -M ^(^) ■N[l]^lf R{z - Uj) ^ 

p1 - , (1) 

where j^Jaj ^ {'-^ i q) denotes the usual nucleon one-body four-current including leading order relativistic contributions 
like Darwin-Foldy and spin-orbit. They contain e.m. Sachs form factors using the dipole parametrization. Explicit 
formulas can be found, for example, in ^Jl- For simplicity, the electric form factor of the neutron is neglected. 
Furthermore, z = W+ie, and = \/M'^ + {p[)'^ + ^/Mjf -f (^2)^ (i^ = 1^1 for any vector a) with W = \/Mj + q^+Lu 
as invariant mass of the photon deuteron system in the cm. frame and Md as deuteron mass. Finally, R denotes the 
so-called dressing operator describing hadronic off-shell effects |^ |^ Q . 

The quantity J^iopsub ™^ O denotes the e.m. corrections arising from the pion-nucleon loop which have been 
discussed in detail in Q and can be interpreted as e.m. off-shell contributions. For the transition to virtual photons 
e.m. form factors have to be incorporated. Moreover, the propagators of the loop diagrams with an intermediate 
7*7riV-state have to be modified due to the fact that w 7^ g in electrodisintegration. 

For the one-body (7-l-iV — > A)-transition current, we take the usual nonrclativistic form including only the dominant 
Ml-transition 

^A\3^N{Wsub,i^,q)\PN) = 5{p/^-Pn -^li) ^^^G^i^{Wsub + i<^,i^,q)iS^N >i q~tN , (2) 

where 

'IjN = q ^VlT^PA with Mr^ = 1232MeV. (3) 

As has been outlined in detail in the e.m. coupling constant has been fixed for real photons by a fit of the 

(3/2) 

M^^ -multipole of pion photoproduction on the nucleon. Due to nonresonant rescattering mechanisms, it becomes 
complex and dependent on the invariant mass Wsub of the TrA'^-subsystem for which we adopt the spectator-on-shell 
approach 3|. Moreover, when embedded into a nuclear medium, G^/( depends also on the photon momentum q and 
the photon energy ui. As has been shown in [12| |. at least for photodisintegration a so-called on-shell prescription 

GfitiWsub + te, u, q ) ^ G^ll{Wsub + ie) (4) 

turns out to be very accurate, simplifying considerably the numerical evaluation. The resulting coupling constant 
used in this paper is parametrized as follows 

G^Hiz = W,ub + - Ji^.N{W,ub) e^*(^=-) (5) 



with 



MAJv(M/snb) = MO + M2 ( — I + M4 — I and -^{Wsub) = ■ r , (6) 



where the on-shell pion momentum q^^ is a function of Wsub according to Wsub ~ + ~^ V + 9? i while for 
Wsub < Mm + m^r we have set = 0. The parameters in ^ have been fitted to the resulting coupling constant 
"G^/^(effl)" discussed in |3] and are listed in Tabled In addition a contribution from the 7A^A-transtion charge 
density is included. For virtual photons, again G^ll must be multiplied with an appropriate e.m. transition form factor 
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which for simpHcity has been chosen in the dipole form. We are aware of the fact that with increasing momentum 
transfer one finds a shghtly stronger fall-off of {Q^) with = —ui^ compared to the dipole form . However, 
we restrict ourselves here to kinematics with < 4 fm~^ where the difference does not matter. 

With respect to the two-body currents, we consider, as in Q, meson retardation in the pure pion exchange con- 
tributions whereas retardation in 77rp/w-MEC turned out to be unimportant in photodisintegration (12j . but very 
CPU-time consuming. Therefore, in order to facilitate our numerical evaluation, we have neglected retardation in the 
latter contribution. Furthermore, we include in addition static p-exchange as well as various A-MEC which have been 
discussed in j^. Similar to the e.m. loop corrections the transition from photo- to electrodisintegration leads to some 
modifications in the intermediate iViV7r7*-propagators as well as for the incorporation of an e.m. form factor. 

In the static model the e.m. loop contributions J'l^^p to the nucleonic one-body current are neglected naturally 

and the hadronic dressing factor R is replaced by 1. Obviously, the MEC contain only the usual static terms. Note 
that in the static limit the recoil contributions do not contribute due to their cancellation against the wave function 
renormalization Moreover, in the static limit there are no nonrelativistic MEC-contributions to the charge 

operator. 

For the numerical evaluation of the T-matrix for deuteron electrodisintegration, we use the standard multipole 
decomposition of the e.m. current. Moreover, we take advantage of Siegert's theorem which allows to express the 
dominant contributions of the transverse electric multipoles via the charge multipoles. 

III. RESULTS AND CONCLUSIONS 

Neglecting polarization effects the differential cross section for deuteron electrodisintegration in the one-photon 
exchange approximation is determined by four structure functions, two diagonal ones /i and /t and two interference 
ones Jlt and /tt They are functions of the squared three momentum transfer in the cm. system, the final 
state cm. -energy Enp = W — 2Mjv, and the angle 9 between q and the relative neutron-proton momentum in the 
final neutron-proton cm. -system. 

We begin the discussion with the experiment of Jordan et al. 16] for the kinematics Enp — 66 MeV, q^ — 3.87 fm^^ 
which has been performed in order to extract the structure functions /l, /t and /lt- Our predictions for these 
structure functions using both static and retarded operators are shown in Fig. ^ together with the few existing 
experimental data points. Since this kinematics is well below pion threshold it is not surprising that both approaches 
yield almost identical results - in fact only for the smallest structure function fxT the corresponding curves can be 
distinguished - demonstrating that both approaches are indeed equivalent in this kinematic region. A similar result 
has also been found in photodisintegration 3]. The agreement with experiment is quite satisfactory but more data 
are certainly needed for a more critical comparison. 

The situation changes completely if the excitation energy is above pion threshold. As a first example, we consider 
in Fig. |21the kinematics Enp = 179 MeV, = 1.66 fm~^ of an experiment by Turck-Chieze et al. [13. Whereas fr 
is only moderately affected by retardation, the other structure functions, in particular and /lt are much more 
sensitive to the inclusion of retardation. At forward angles /l is enhanced by roughly 45 percent while above 6 « 90° 
it is significantly lowered. Especially Jlt increases strongly in absolute magnitude. A detailed analysis has shown 
that this strong retardation effect is mainly due to the charge recoil contribution to the Coulomb monopole (see Fig.|2| 
for a graphical representation). This is demonstrated in Fig.|21by the dash-dotted curves for which this contribution is 
switched off. The remaining retardation effect is significantly smaller. A comparison with the experimental diffferential 
coincidence cross section is shown in Fig. 0] One readily notes a considerable enhancement of the cross section in the 
retarded approach improving significantly the agreement with the data. But the agreement is not perfect, the theory 
being a little too low in the maximum around 9^^^ = 45° and somewhat too large at higher angles. Again, the most 
important contribution of retardation is due to the recoil charge contribution to the monopole. 

The importance of retardation becomes even more striking for excitation energies Enp in the A-region. This is 
demonstrated in Figs. [S] and for the kinematics of a more recent experiment by Pellegrino et al. with Enp = 280 
MeV and — 2.47 fm^^. The effects of retardation on anf f^r in Fig.|Slare even stronger than in the previous 
example. However, in the differential cross section the retardation effects in the different structure functions appear 
to cancel each other to a large extent so that static and retarded approaches yield very close results in the measured 
angular range 95° < 6 < 135° (see Fig. ^ but lead to a slight underestimation of the experiment. The better 
agreement between theory and experiment reported in 18] is obtained by the use of a modified 7iVA-coupling which 
is considerably stronger than ours given in Table ^ due to the assumption of a vanishing nonresonant contribution 
to the M{i^^^''-multipole. However, as has been discussed in this approach causes formal inconsistencies and is, 
therefore, questionable. 

Comparing the static as well the retarded results for /tt in Fig. to experimental data, one notes a sizable 
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underestimation of both approaches. This discrepancy may have different origins. First of all, one has to be aware of 
the fact that the majority of the data (full circles in the right panel of Fig.|SJ) has been extracted under the assumption 
Jlt — which is well justified only in the static, but not all in the retarded approach. Only the open square data 
points have been obtained without this assumption. Concerning the theoretical uncertainties, the present approach 
does not contain all relativistic contributions of leading order in p/M^ as is the case in the approach, for example, of 
Ritz et al. jl9l |. On the other hand, among the relativistic contributions retardation effects are unique with respect 
to the singularity structure of the corresponding meson-nucleon propagators above pion threshold which cannot be 
treated within a Taylor expansion in p/Mj^. Therefore, concerning the role of retardation, the approach of (IQj cannot 
be considered realistic above pion threshold. However, the model of Ritz et al. is useful to estimate the role of boost 
effects (which do not contain any singularities) which are still missing in the present approach. It turns out that for 
the kinematics of Fig. El they are small. Moreover, we have checked that additional, not yet included mechanisms like 
static (T- or ^-MEC, a nonvanishing Geu or a possible £'2-excitation of the A are small and cannot explain the noted 
discrepancy between theory and experiment in Fig. |B1 

In summary, we have studied the role of retardation effects in electrodisintcgration of the deuteron for various 
kinematics. It turns out that especially the recoil charge contribution to the monopole, which is not present in 
static approaches, is very important for excitation energies above pion threshold, leading to dramatic changes in the 
structure functions /l and Jlt whereas the other structure functions fx and fxT are much less affected. 

With respect to future developments, one important step concerns an improved treatment of the final-state- 
interaction. As has already been noted in 0], the description of the phase shifts in several NN-partial waves is 
only fairly well above pion threshold for our NN-interaction model. This touches a general problem of realistic NN- 
interactions: In contrast to the energy regime below pion threshold, at present no "high precision" NN-interaction is 
available for energies above pion threshold. However, this is an indispensable ingredient for any theoretical description 
of e.m. reactions on the deuteron above pion threshold like photo- or electrodisintcgration, but also meson production 
on the deuteron. A further interesting topic will be the study of polarization observables which supposedly are more 
sensitive to interaction effects. 
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TABLE I: Parameters of the A-coupling in 0. 
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FIG. 1: The structure functions Jl, fr, Jlt and Jtt for Enp = 66 MeV, q 
static approach; full: retarded approach. Experimental data from fl^ . 
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FIG. 2: The structure functions /t, /lt and Jtt for i5„p 
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The additional dash-dotted curves represent the results of the retarded approach where the Coulomb monopole contribution of 
the recoil charge operator is switched off. 




FIG. 3: Diagrammatic representation of the recoil contributions to the effective charge and current operators. The coupling 
of the photon to the nucleon is given solely by the nonrelativistic charge and spin/convection current, respectively. 




FIG. 6: The differential cross section (left panel) and the /TT-structure ftmction (right panel) for E„p = 280 MeV, 
2.47 fm"^,6'e = 30°. Notation as in Fig.[Tl Experimental data from (Q. 



